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Flow signatures in experimental data from relativistic ion collisions are usually interpreted as a 
fingerprint of the presence of a hydrodynamic phase during the evolution of these systems. In this 
work, flow signatures arising from event-by-event viscous hydrodynamics are compared to those 
arising from event-by-event non-interacting particle dynamics (free-streaming), both followed by a 
late-stage hadronic cascade, in d+Au, 3 He+Auat yfs = 200 GeV and p+Pb collisions at y'i = 5 
TeV, respectively. For comparison, also Pb+Pb collisions at yfs = 2.76 TeV are simulated. It is found 
that non-hydrodynamic evolution can give rise to equal or larger radial flow than hydrodynamics 
with rj/s = 0.08 in all simulated collision systems. In light-on-heavy-ion collisions, free-streaming 
gives rise to triangular and quadrupolar flow comparable to or larger than that from hydrodynamics, 
but it generally leads to considerably smaller elliptic flow. As expected, free-streaming leads to 
considerably less elliptic, triangular and quadrupolar flow than hydrodynamics in nucleus-nucleus 
collisions, such as event-by-event Pb+Pb collisions at y'5 = 2.76 TeV. 
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I. INTRODUCTION 


One of the recent great successes of high-energy nuclear physics has been the understanding that the hot QCD 
matter created in high-energy nucleus-nucleus collisions behaves like a strongly-coupled, low-viscosity fluid rather 
than a weakly-coupled, almost ideal gas of quarks and gluons. This finding rests on the fact that several, independent 
observational probes of the hot QCD matter, such as for example the strong collective flow and strong jet quenching, 
can all naturally be explained theoretically under the single assumption of a hydrodynamic phase early in the evolution 
of the hot QCD matter. 

However, the sizable flow signals that have been measured in systems created in light-on-heavy-ion collisions at 
RHIC and the LHC [THS] have come as a surprise to many experts in the field. There are currently two competing 
interpretation for this finding. On the one hand, the flow signals measured in light-on-heavy-ion collisions may 
just have the same origin than those in nucleus-nucleus collisions, namely the presence of a hydrodynamic phase 
(see Refs. [7HT4] for work along those lines). Because the systems created in light-on-heavy-ion collisions are very 
small and short-lived compared to those created in nucleus-nucleus collisions, one may, however, entertain a different 
hypothesis, namely that flow is being created by some other, non-hydrodynamic mechanism (see Refs. [T5HT5] for 
work along those lines). 

In view of these competing interpretations, the present work tries to answer the following question: How different 
would flow observables be if the systems created in relativistic ion collisions never go through a hydrodynamic phase, 
but still experience expansion and cooling in the hot phase and interactions in the low-temperature, hadron gas phase? 

To answer this question, the dynamics of the hot QCD matter phase is alternatively described by two extreme 
(classical) opposites: hydrodynamics and a non-interacting gas of free-streaming particles. By keeping all steps of the 
model simulations the same but only switching from a hydrodynamic to a free-streaming description (with the same 
equation of state) allows one to make fully transparent comparisons between hydrodynamic and non-hydrodynamic 
results. The extreme nature of the two opposites employed also guarantee that any intermediate case between strong 
and weak interactions must be bounded by the results found in this study. 

It is important to stress the relation of the present study to previous works. For instance, Refs. mmsi proposed an 
initial state effect originating in QCD to explain the observed elliptic flow in light-on-heavy-ion collisions, a mechanism 
that seems to be in tension with more recent experimental results for multi-particle correlations 0. This line of work 
is complementary to the present study, since in the present study the emphasis is on non-hydrodynamic flow generated 
during the hot QCD phase 

Refs. [TTl jT8] describe the hot QCD matter phase using the phenomenological, non-hydrodynamic AMPT model 
including strings, hard particles and effective interactions, as well as hydrodynamics. AMPT seems to match the 
experimental results from p+Pb collisions at y?s = 5.02 TeV with a very small effective scattering cross-section, 
thus suggesting the possibility a non-hydrodynamic explanation for the observed flow signal. However, the fact that 
AMPT has many different phenomenological ingredients makes the interpretation of this result somewhat challenging. 
The main difference to Refs. [ED HE] of the present study is the fully transparent nature of the non-hydrodynamic 
description employed here. 

Finally, Ref. [T9] uses the hadronic cascade model URQMD to describe p+Pb collisions at y/s = 5.02 TeV. It was 
found that purely hadronic interactions cannot describe experimental data in p+Pb collisions, but no comparison 
with hydrodynamics using the same equation of state as in the URQMD model was attempted. The present study 
is similar to Ref. m in that a hadronic cascade code is employed to describe the low temperature system evolution. 
A key differences to Ref. [19] in the present study is that besides purely hadronic evolution (applicable to the low 
temperature phase), also two different scenarios for the hot QCD phase (hydrodynamics and free streaming) are 
considered. 


The remainder of this work is organized as follows: In Section [IT] the details for the model setup, initial conditions, 
equation of state, and calculation of final observables are given. Also, this section contains a warm-up example of 
comparing free-streaming and hydrodynamics in the case of collisions of smooth nuclei. In Section m results from 
simulating event-by-event collisions of granular light nuclei on heavy nuclei (p+Pb, 3 He+Auand d+Au), as well 
as nucleus-nucleus collisions (Pb+Pb) are shown. The summary and conclusions are presented in Sec IV, while the 


special case of an analytically solvable comparison between ideal fluid dynamics and free-streaming is discussed in 
appendix [X] 
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II. METHODOLOGY 

A. Stage 1: Initial Conditions and Equation of State 


Initial conditions are prepared for the bulk energy-momentum tensor T ab . For simplicity, only boost-invariant 
dynamics will be considered, which is most easily implemented by using an expanding space-time described through 
the Milne coordinates: 


T = \Jt 2 - Z 1 , 



t + z 
t — z ’ 


X_L 



(1) 


For a boost-invariant system, the dynamics is invariant under translations in £, so that all quantities only depend on 
t, Xj_5 thus, initial conditions amount to specifying T ab {r = to,xj_) at some time tq. 

Semi-realistic initial conditions for T ab have been calculated under certain assumptions of the collision dynamics 
[2(M(3| . The collection of these amount to a large class of (generally non-equilibrium) initial conditions that could be 
implemented for T ab (r = To). The goal of this study, however, is the comparison between the subsequent evolution 
of T ab (r, xj_), and thus arguably the simplest possible initial condition will be implemented in the following: the case 
of thermal equilibrium with zero local flow in the transverse x_l plane. In this case, the initial conditions are fully 
specified by the energy density T 00 (xj_) and the equation of state (e.g. through the functional relation of the pressure 
P(e) to the energy density e). It should be stressed that an equilibrium initial condition (or more precisely an initial 
condition with zero momentum anisotropy) is not required for the comparison. Also, it should be pointed out that 
a condition with zero momentum anisotropy is generated ’by accident’ at some point during the time evolution for 
initially prolate momentum distributions (cf. 1271 1. In this sense, the initial conditions chosen above do not actually 
require equilibration of the system to happen before r = To, but correspond to a whole class of non-equilibrium initial 
conditions at specially chosen instances in time. 

Within hydrodynamics, arbitrary equations of state are easily implemented. In the case of classical particle dynamics 
described through a Boltzmann equation, this is considerably harder to do (see e.g. Ref. [28] on how to implement 
generic equations of state through generalized Boltzmann equations). However, since the present study aims at 
mapping the hot QCD dynamics onto a hadron gas cascade evolution at low temperatures, the absolute minimum 
requirement in order to conserve energy and momentum throughout the whole evolution is to implement an equation 
of state that smoothly matches onto the equation of state from a hadron gas at some predefined switching temperature 
T = Tgw ■ This can be achieved by considering the equation of state generated by a number of Z particle with mass 
to: 


e 


Z m 2 T 
27r 2 


(»■*,(=) + m*(=)) 


P = 


Z m 2 T 2 
2n 2 


K 2 



( 2 ) 


Choosing to match at Tsw = 0.17 GeV, it is found that the hadron gas pressure as well as its first derivative can be 
matched at T = Tsw by choosing to = 0.779 GeV and Z = 116. Figure [l] shows a comparison for the pressure from 
the massive gas, the pressure from the sum of all known hadron resonances up to masses of 2.2 GeV in the Particle 
Data Book, and the pressure calculated in lattice QCD by the BMW collaboration [25]. As can be seen, the fact 
that the value of the pressure as well as its first derivative match between the hadron gas equation of state and the 
one-component gas imply that a smooth transition from one description to the other is possible. For this reason, the 
one-component gas equation of state with parameters m = 0.779 GeV and Z = 116 will be adopted for the rest of 
this work, both in the non-interacting and hydrodynamic case. 

One could worry that in the case of free-streaming the system would generically be far from equilibrium, and an 
equation of state for such as system could not be defined. However, non-perturbative quantum field theory studies 
provide evidence that an equilibrium equation of state relating the energy density to the temperature is present 
after a very short time-scale even for otherwise out-of-equilibrium situations [3(31. Thus, even far from equilibrium 
systems can reasonably be expected to possess a relation between energy density and temperature that is given by 
the equilibrium equation of state. 

The initial conditions for the energy density are taken from two sources: first, a smooth optical Glauber model 
which is used to demonstrate key features of the free-streaming results in a simple setting. And second, from a more 
realistic event-by-event Monte-Carlo Glauber event generator taking into account the collision system composition 
and collision energy. In the first case, the different nuclei are modeled by employing an overlap function 


dz 


1 _|_ g-(l x i| 2 +i 2 -fl)/a 


Ta(x, y) 


(3) 
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Equation of State: Pressure 



FIG. 1. Pressure versus temperature for different equations of state. Shown are results for the hadron resonance gas and the 
one-component massive gas with parameters adjusted to match the value and slope of the hadron gas pressure at T = Tsw- 
Inset shows dependence of speed of sound square c 2 . For illustrative purposes, the full lattice QCD result (from Ref. [29]) is 
also shown. 


with R, a the charge radius and skin depth parameters (in the following, R = 4.16 fm and a = 0.606 fm for 63 Cu will 
be used). The initial condition for the energy density in this case is then given by 

e0b,x_L) =E q T a (x + ^,y)T A (x - ^,y), (4) 

where b controls the impact parameter of the collision and Eq is an overall constant controlling the total multiplicity 
(entropy) simulated. 

In the second case, initial conditions for the energy-density e(ro,xj_) for each event are constructed as follows. 
Using Woods-Saxon distribution functions for the heavy ions such as Au, Pb mm, the Hulthen wavefunction for 
the deuteron (cf. [32]) and realistic calculations for the 3 He wavefunction [33], probability distributions of the nucleons 
within the nuclei of interest (cf. ]Sj) are obtained. Using a Monte-Carlo Glauber [35], these probability distributions 
are mapped to positions of individual nucleons in the transverse (x, y) plane on an event-by-event basis implementing 
a hard-core repulsive potential of radius 0.4 fm between nucleons. The positions of nucleons undergoing at least one 
inelastic collisions are recorded (“participants”) and converted into a density function i? 2 (x) by assuming that each 
participant contributes equally as a Gaussian with a width of w = 0.4 fm (to match the RMS radius of a single 
nucleon). The initial condition for the energy density is then assumed to be given as 

e ( r 0j x ) = E 0 R 2 (x) , (5) 

with Eo again an overall constant (dependent on to, collision energy and collision system) that is related to the total 
multiplicity of the event. Typically 100 initial conditions are generated for each collision system. 


B. Stage 2, Option a: Hydrodynamics 

Once the initial conditions for the energy-momentum tensor are specified, these can be converted into hydrodynamic 
degrees of freedom via the following decomposition of the energy-momentum tensor: 

T ab = eu a u b - (P - n) (g ab - u a u b ) + n ab , (6) 

where e, P are the local (equilibrium) energy density and pressure, u a is the local fluid four velocity, g ab is the metric 
tensor and II, ir ab are the shear and bulk stress tensors, respectively. Such a decomposition of the energy-momentum 
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tensor is possible in most cases, with the exception of out-of-equilibrium quantum states (see e.g. Ref [551 for more 
discussions). For the case of equilibrium conditions at hand, the initial conditions imply u a = (1,0, 0, 0), II = 7 r° & = 0 
and e(T 0 ,XjJ given by the initial conditions discussed in Sec IIA 


Once the initial conditions have been specified, the hydrodynamic equations of motion V a T ab = 0 have to be solved. 
To do this, one first needs to specify the constitutive relations that e.g. connect the shear tensors to gradients of 
the fundamental hydrodynamic degrees of freedom e, u a . Fortunately, recent progress in relativistic fluid dynamics 
(which to a large extent has been fueled by access to strongly coupled field theory dynamics from the gauge/gravity 
duality conjecture) has led to a complete characterization of all possible terms that can appear to a certain order in 
gradients (see e.g. Refs. [37H40] ). Once the constitutive relations have been specified, one still needs an algorithm to 
actually solve the hydrodynamic equations of motion numerically. The standard approach in relativistic dissipative 
fluid dynamics, which has been fully developed in the past 10 years, is to use causal second-order fluid dynamics (see 
Refs. pl»45] for reviews on this subject). 

Finally, there are transport coefficient functions appearing in the hydrodynamic equations of motion. At zeroth 
order in gradients (ideal fluid dynamics), the only such quantity is the speed of sound c s , which is fully specified 
through the equation of state c 2 = At first order in gradients (Navier-Stokes fluid dynamics), there are the shear 
and bulk viscosity coefficients, denoted as 77 , £, respectively. For simplicity, in this study only constant values for the 
ratio of shear viscosity over entropy density s will be adopted, and the bulk viscosity will be set to zero. Finally, at 
second-order in gradients there are an additional 11 transport coefficient in flat space-times |59j . and for simplicity 
most of these will again be set to zero, except for the relaxation time r n = 

With these specifications, the equations of motion are solved using the publicly available VH2+1 code package 
.Us TTj . version 2.0, on a two-dimensional space grid with lattice spacing of Ax ~ 0.1 fm. Every 0.25 fm/c during the 
evolution, the local temperature of fluid cells is monitored and once a fluid cell cools below the switching temperature 
Tsw ; information about the cell’s location as well as the value of e, u a , II, ir ab is stored. The collection of all these cells’ 
locations (t, xj_) defines the switching hypersurface X, which will eventually be used to initialize the low-temperature 
hadron gas dynamics (see IID). 


C. Stage 2, Option b: Free-Streaming 


While the above hydrodynamic option to describe the bulk system dynamics is quite standard, this work proposes 
an “option b” for the dynamics: non-interacting free particle dynamics. In this case, the energy-momentum tensor is 
given in terms of the one-component on-shell particle distribution function /(r, x_l, £, pj_,p^) as 


T 


ab 


d 2 p±dp^r p a p b 
(27t) 3 p T 


f(r, xj_,£,p_l,/), 


(7) 


where for on-shell massive particles p T = m 2 + p\ + r 2 p f 2 . The distribution function will be taken to be a solution 
to the classical Boltzmann equation in the non-interacting (free-streaming limit) [281 : 

p a d a /-—d' p) / = 0, (8) 

T s 

where d^ This equation is readily solved using the method of characteristics, finding the general solution 


/ = / 


^Pj_,P 5 ,x_l 


rp±p T 

2 1 9 ’ 

p + m z 


£ + ln 


V 

.pc 



(9) 


Also, the implementation of equilibrium initial conditions is straightforward. Given an equation of state, the energy 
density defines a local equilibrium temperature T = T(e) and an equilibrium solution for the particle distribution 
function for Eq. Q can be shown to be given by 


/ eq = Ze~P au ^ T , 


( 10 ) 


where u a is the local macroscopic (not necessarily fluid) four velocity with respect to some global laboratory frame, 
and Z is the effective number of degrees of freedom first introduced in Eq. §■ It is straightforward to show that 
inserting Eq. (10 1 into Eq. 0 leads to the results given in Eq. 0. Evaluating Eq. ( [To| for T = T (e(xj_)) in the 
transverse plane then fully specifies the initial conditions for the free-streaming dynamics. 
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For the case of boost-invariant dynamics, and equilibrium initial conditions with m°(to) = (1,0, 0, 0) given at r = To, 
the solution to Eq. |8]) at any later time may then be analytically written asQ 


/(t>x _l,£, P±,/) = Z exp 


Po/ T ( x ^ - P± ^[+m^ ) ’ ^ = ^P‘i + m2 +pV t $ ■ ( n ) 


From the solution at time r, one can evaluate the energy-momentum tensor T ab and from the energy-momentum 
tensor one can find the local energy density, flow velocity, shear and bulk stress tensors using the decomposition 
in Eq. ©• (Note that since the particle dynamics is classical, a decomposition along the lines of Eq. (Jb]) is always 
possible even for far-from-equilibrium systems [31]). Using the same routines as in the hydrodynamic framework, the 
local temperature is monitored and a switching hypersurface E can again be defined as those space-time points which 
have T = Tsw (in practice, and for better comparability, the same two-dimensional lattice as in the hydrodynamic 
framework with Ax ~ 0.1 fm and a time-increment between steps of 0.25 fm/c is used). The quantities e, u a , II, n ab are 
stored along the hypersurface and thus the final information available is exactly equal to that from the hydrodynamic 
framework. 


D. Stage 3: Kinetic Freeze-Out and Hadron Cascade 


Using information from the switching hyper-surface from either the hydrodynamic or free-streaming evolutions 
in the hot phase, the low temperature phase is simulated through a hadronic cascade code (B3D, [55]). Using the 
hyper-surface information to boost to the rest frame of each cell, the cascade is initialized with particles in the rest 
frame drawn from a Boltzmann distribution at a temperature Tsw with modifications of the momentum distribution 
to include deformations from viscous (both shear and bulk) stress tensors (see [55] for details). Specifically, for a 
particle with mass [i, the distribution in the local rest frame is assumed to be of the form 


/(P) = exp 


-vV+k 2 /T' 


Pi — (&ij “b Ajj ) Pj i 


( 12 ) 


with A ij controlling the size of the shear and/or bulk corrections to the stress-energy tensor. Note that for A tJ <C 1, 


one finds /(p) — /eq(P) oc f eq (p) 


plpjX ' J [46] with T’ = T. However, even for small values of the stress tensors, the 

because of the additional power of p^u^ in the 


/eqUV ^ Je q vr , pQ 

distribution function does not recover ’quadratic ansatz’ form 
denominator resulting from expanding Eq. (12). It should be stressed that this procedure ensures that the complete 
stress tensor T ab (not just its ideal fluid part) is matched across the hypersurface boundary. In particular, this implies 
that no assumption about equilibrium is made at the switching hypersurface: arbitrary deviation from equilibrium, 
parametrized through large dissipative tensor components, are allowed. If only the ideal fluid part of the stress-energy 
tensor is matched, this would lead to a ’fake’ collective flow signal, such as a discontinuous elliptic flow component 
(cf. [57]). See the discussion in Appendix |B| for more details. 

The cascade code B3D includes hadron resonances in the Particle Data Book up to masses of 2.2 GeV, which 
interact via simple s-wave scattering with a constant cross-section of 10 mb as well as scattering through resonances 
(modeled as a Breit-Wigner form). Once the resonances have stopped interacting, one can obtain final charged hadron 
multiplicities d ^y h , mean charged particle momentum (j?t) and flow coefficients v„ (pt) for n > 1 from summing over 
individual particles with momenta p. Specifically, 


dN ch 


E 


,ch. particles 
PT 


2 n pTdYdpT 2 irpT At Ay 


dNch. 

dY 


= fo°°dP T 


dN c h 
dY dpT 


(pr) = ^ 


dp T PT 


dN cL 

dY 


dY dp'j' 


K|(Pt) = \/Snipr ) 2 + C n (p T ) 2 , 


Sn{pr) \ = 
Cn(pr) ) 


E 


sin (n<j))) 
cos (ncj)) 

v >ch) particles 


^ch. particles 
PT bin 


<p = arctari ( T ^) , (13) 


where A t = 80 MeV, Ay = 2 are the width of bins for particle pr and rapidity Y, respectively. Note that since 
the cascade is applied to a boost-invariance case, the large Ay value is of no significance. In practice, a sum over 
both particles and anti-particles and division of the spectra by two is performed, in order to increase statistics. For 
every hydrodynamic evolution event, at least 100,000 B3D events are run to increase statistics. In doing so, the 
sums in the definition of v n above are extended over all B3D events, thereby explicitly ignoring fluctuations arising 


1 Note that the result correspond the form used in anisotropic hydrodynamics m ■ 
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from hadronic decays. After thus obtaining results for 2 TT PT dYdp T an< ^ v n(PT ) for each hydrodynamic event, an event 
average to obtain the event-by-event mean and event-by-event fluctuation is performed. Results both for the case 
of hydrodynamic and free-streaming hot phase dynamics are reported on in the following. It should be noted that 
because B3D enforces detailed balance, there is no baryon annihilation simulated and as a consequence proton yields 
are too high. For this reason, the results reported for protons below should be interpreted with care. 


E. Warm-up: Collisions of idealized smooth nuclei at b = 4 fm 

As a warm-up example, consider the case of smooth optical Glauber initial conditions for collisions of 63 Cu nuclei 
at an impact parameter of b = 4 fm. In this case the initial conditions are simple enough that the main physics 
similarities and differences between hydro and non-interacting gas can be understood. 

In Fig. [2] time-snapshots of the temperature, velocity and shear tensor space profiles along the y-axis are shown. 
One notes that despite the very different character of the hydrodynamic and free-streaming evolution, the temperature 
profiles during most of the evolution are almost identical. The equal-time velocity comparison shows that flows are 
also similar in magnitude, but the velocities from non-interacting evolution are consistently larger than those from 
hydrodynamics. This is easy to explain: in almost ideal hydrodynamics, the pressure along the transverse axes 
Pt and the longitudinal axis Pl are almost identical (after all, hydrodynamics implies that the system is locally 
approximately isotropic). By contrast, in the free-streaming evolution in the boost-invariant approximation the 
longitudinal pressure falls quicker than the transverse pressure because there are no particle interactions to keep the 
system locally isotropic. Since the sum of the transverse and longitudinal pressure is fixed by the equation of state, 
this implies that the transverse pressure in free-streaming will generally be larger than the transverse pressure in 
almost ideal hydrodynamics. Since flow velocities in the transverse plane are being sourced by the gradient of the 
pressure, larger transverse pressures lead to larger flow velocities in the non-interacting case. For a particular case 
where the similarity of radial flow in free-streaming dynamics and hydrodynamics can be analytically demonstrated 
see the discussion in Appendix |A| 

When considering the result for the shear stress tensor in Fig. [2] the above similarities between hydrodynamic and 
free-streaming evolution stop. While the shear stress is generated gradually in hydrodynamics, the non-interacting 
free-streaming evolution leads to a sudden build-up and subsequent saturation of the shear stress. However, it is 
interesting to note that despite the fact that the free-streaming evolution literally corresponds to infinite viscosity, the 
overall magnitude of the shear stress tensor generated during the evolution is comparable to that from hydrodynamics 
with extremely small viscosity over entropy ratio 77 /s ~ 0.08. 

The information about the components of the energy-momentum tensor is imprinted onto the final particle spectra, 
shown also in Fig. [2j It is important to recall that the final particle spectra result from either hydrodynamics or 
free-streaming dynamics in the hot phase of the evolution T > T$w followed by the same hadronic cascade evolution 
in the cold phase T < Tsw- Shown are results for pions, kaons and protons, and the larger transverse flow developed 
in the free-streaming evolution (as compared to hydrodynamics) is clearly seen as a flattening of the spectra for all 
particle species. From this figure, it is evident that radial flow does not indicate the presence of a hydrodynamic phase 
during the system evolution. This has been noticed before [22] • Nevertheless, it is worth stressing that the presence 
of radial flow should not be used as an indicator for hydrodynamics, as has been often assumed (see e.g. Ref. [H]). 

Also shown in Figj2] is the identified particle elliptic flow (pr) resulting from hydrodynamics or free-streaming 
dynamics. Note that in the simple case of smooth initial conditions, by symmetry this is the only non-trivial anisotropic 
flow component v n . From this plot it is evident that almost ideal hydrodynamics gives rise to a considerably larger 
elliptic flow than free-streaming dynamics, confirming the predominant view that hydrodynamics is necessary to 
explain strong anisotropic flow. However, the elliptic flow found for the case of free-streaming is not consistent with 
zero. At first glance, this is puzzling, given that it can be analytically shown that free-streaming does not generate 
momentum anisotropies by itself, while diluting the spatial anisotropies (see e.g. Refs. [441 [49]). However, even 
though the spatial anisotropies are being diluted, their potential to generate momentum anisotropies is not actually 
lost. Rather, what is happening is that both macroscopic velocities and dissipative parts of the stress tensor are 
being generate in precisely such a way that the net (non-equilibrium) momentum anisotropy is exactly zero. This is 
demonstrated explicitly in AppendixjB] In the case at hand, the full energy-momentum tensor after the free-streaming 
evolution is used to initialize the late-stage hadronic evolution, and it turns out that the hadronic interactions are 
sufficient to re-generated part of the momentum anisotropies from this T ab by strongly damping the dissipative parts 
while the flow velocities remain. This is only possible if the hadronic evolution itself has transport properties similar 
to a “low” viscosity fluid, because otherwise momentum anisotropies on the level seen in Fig. [2] (e.g. 50 percent of 
hydrodynamics with 77 /s = 0.08) would never be (re-)generated. Recent measurements of the ratio of shear viscosity 
over entropy density in the hadron gas phase are consistent with this picture m- 

The presence of a hadron gas phase (often referred to as “corona” in earlier work |511 1 is essential for generating 
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FIG. 2. Free-streaming evolution (no-interaction) and almost ideal hydrodynamics (r)/s = 0.08), followed by a hadronic 
cascade for smooth Cu+Cucollisions at b = 4 fm. Shown are time-snapshots of the temperature profile (upper left), velocity 
profile (upper right), shear tensor profile (middle left), switching hypersurface (middle right), identified particle spectra (lower 
left) and identified particle elliptic flow (lower right). For reference, the switching temperature is indicated in the upper left 
plot. Note that the ’true’ theoretical V 2 (pr) curves would be smooth, but the finite sampling statistics for the hadron cascade 
code introduces some statistical error that is particularly evident for the smaller values encountered in the free-streaming case. 
See text for details. 


the anisotropic flow effects seen in the identified particle plots in Fig (2j Without hadron gas phase, there would be 
radial flow (see the analytic result presented in Appendix [A}, but no elliptic flow. However, in actual systems created 
in relativistic ion collisions there always is a hadronic gas phase, so it is crucial that this component be included in 
the system description, and that its transport properties are better quantified (see e.g. Ref. [SO] for work along these 
lines). 
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Particle Spectra: Hydro versus Non-Interacting Gas 
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FIG. 3. Simulations of granular p+Pb collisions at yfs = 5 TeV. Shown are final particle spectra and anisotropic flow coefficients 
Vn (pt) for identified particles for free-streaming evolution (no-interaction) and almost ideal hydrodynamics (r)/s = 0.08), 
followed by a hadronic cascade. See text for details. 


III. RESULTS 


A. Central event-by-event collisions of granular nuclei 


A more realistic application of the techniques highlighted in the previous section is the relativistic collision of light- 
on-heavy-ions, such as 3 He+Auand d+Auat ^J~s = 200 GeV and p+Pb at y/s = 5.02 TeV. For each of these collision 
systems, one hundred initial events are generated from a probability distribution of nucleons inside the colliding nuclei 
(see section IIA for a detailed discussion). For each of these events, the subsequent dynamics is simulated using either 


a hydrodynamic evolution or a free-streaming evolution, followed by the same hadron cascade for the low temperature 


phase. Unlike the simplified case discussed in the section HE the granular nature of each individual event gives rise 
to all anisotropic flow harmonics v n with n > 1, not just the elliptic flow V 2 - 

The results for p+Pb collisions at a/s = 5.02 TeV are shown in Fig. [ 3 ] Considering the identified particle spectra, one 
finds that the additional radial flow generated in the free-streaming dynamics compared to hydrodynamics is almost 
negligible, and the resulting spectra are essentially indistinguishable. One reason for this may be the comparatively 
shorter evolution time spent in the hot phase T > T$w for p+Pb collisions compared to the case of smooth nucleus- 
nucleus collisions considered in Sec III El 

The comparison between free-streaming dynamics and hydrodynamics for the elliptic flow coefficient V 2 are con¬ 
sistent with the findings for smooth nucleus-nucleus collisions considered above: the coupled free-streaming and 
hadron gas dynamics gives rise to a non-negligible amount of V 2 , but it is considerably less than the V 2 generated in 
hydrodynamics. 

Considering the higher flow harmonics t> 3 ,V 4 , the comparison between free-streaming and hydrodynamics reveals 
that it becomes more difficult to distinguish between the two scenarios in terms of flow magnitude. For instance, the 
i >3 found for free-streaming plus hadron cascade dynamics is very similar in magnitude to that for hydrodynamics 
plus hadron cascade. Maybe even more interesting, the V 4 amplitude for the free-streaming plus cascade simulation 
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FIG. 4. Simulations of granular d+Au collisions at yfs = 200 GeV. Shown are final particle spectra and anisotropic flow 
coefficients v n (pr) for identified particles for free-streaming evolution (no-interaction) and almost ideal hydrodynamics ( r//s = 
0.08), followed by a hadronic cascade. See text for details. 


in p+Pb collisions turns out to be larger than the corresponding result from hydrodynamics with g/s = 0.08 (see 

Fig.§. 

Overall one finds that free-streaming dynamics followed by hadron cascade dynamics generates approximately the 
same magnitude of anisotropic flow for i> 2 , V 3 and V 4 , e.g. independent from the order of the harmonic. This is clearly 
very different from hydrodynamics, where successively higher orders are more strongly suppressed. 

The results obtained for the p+Pb collisions at y/s = 5.02 TeY should be compared to the results for d+Au and 
3 He+Au collisions at = 200 GeV energies shown in Figs. 4|5 Overall, the same trends that were identified in 
p+Pb collisions repeat for these lower-energy collision systems. However, in d+Au and 3 He+Au collisions at yfs = 200 
GeV one finds that free-streaming plus hadron cascade dynamics generates larger 113,114 than hydrodynamics with 
r)/s = 0.08. Only final V 2 is larger in hydrodynamics than in free-streaming. 

Finally, it is curious to note that the proton 113 is much smaller than pion and kaon V 3 in free-streaming plus 
cascade dynamics in p+Pb , d+Au and 3 He+Au collisions compared to the case of hydrodynamics plus cascade. This 
could suggest a potential experimental handle on separating 113 generated by hydrodynamics from 113 generated by 
non-hydrodynamic processes such as free-streaming. It should be cautioned that simulated proton results could be 
unreliable because of the missing baryon annihilation process in the hadron cascade, see SeclHDl Nevertheless, the 
fact that a large proton V 3 suppression is seen in free-streaming plus cascade but not hydrodynamics plus cascade 
using the same cascade code and for almost identical final identical final particle spectra seems to suggest that this 
effect could be robust. 


B. Mid-central event-by-event nucleus-nucleus collisions 


In view of the above findings for relativistic light-on-heavy-ion collisions, a comparison to nucleus-nucleus collisions 
is in order to see if similar effects are found in these larger systems. Specifically, Pb+Pb collisions at y/s = 2.76 TeV 
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FIG. 5. Simulations of granular 3 He+Au collisions at y/s = 200 GeV. Shown are final particle spectra and anisotropic 
flow coefficients v„(pt) for identified particles for free-streaming evolution (no-interaction) and almost ideal hydrodynamics 
( r//s = 0.08), followed by a hadronic cascade. See text for details. 


and 30%-40% centrality are studied. Results are shown in Fig. [dj The particle spectra are broadly consistent with 
the warm-up case studied in section |II E[ indicating a slightly larger radial flow in the free-streaming dynamics than 
in hydrodynamics. The elliptic flow coefficient generated in free-streaming dynamics of Pb+Pb collisions is found to 
be much smaller than the elliptic flow generated in hydrodynamics at - = 0.08. For higher flow harmonics ( 113 ,^ 4 ) 
the ratio between hydrodynamics and free-streaming is even larger than for V 2 ■ This strongly indicates that for 
larger systems, the longer time spent in the hot (T > Tsw) phase prevents the free-streaming dynamics to preserve 
the space-anisotropies until the hadron cascade takes over, effectively leading to a strong decrease of momentum 
anisotropies with respect to almost ideal hydrodynamics. In essence, this confirms the established paradigm that the 
magnitude of anisotropic flow measured in nucleus-nucleus collisions requires a hydrodynamic phase be present during 
the system evolution. 


One could be worried that this conclusion could be avoided by shortening the time spent in the free-streaming phase 
through increasing the switching temperature Tsw- However, note that for this study Tgw = 170 MeV was chosen. 
Increasing Tsw even further seems to not be justifiable since a hadron gas description is disfavored from lattice QCD 
calculations of the equation of state for T > 170 MeV (cf. Ref. [25]). 


Also, because the anisotropic flow signals in hydrodynamics dwarf the free-streaming results for mid-central 
Pb+Pb collisions at yfs = 2.76 TeV in Fig. [ 6 ] one expects qualitatively similar results for mid-central Au+Au collisions 
at y/s = 200 GeV. However, either in very low collision energy or very peripheral nucleus-nucleus collisions the time 
the system spends in the hot QCD phase is presumably comparable to that in central light-on-heavy-ion collisions, 
so that for these systems one can expect results along the lines of FigsIM 
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FIG. 6. Simulations of granular Pb+Pb collisions at y/s = 2.76 TeV. Shown are final particle spectra and anisotropic 
flow coefficients v„(pt) for identified particles for free-streaming evolution (no-interaction) and almost ideal hydrodynamics 
( r//s = 0.08), followed by a hadronic cascade. See text for details. 


C. Pion Femtoscopy 


Besides flow signals, other experimentally accessible signals such as femtoscopic measurements are often used to 
infer the presence of a hydrodynamic phase in the evolution. 

In this work, the femtoscopic measurements are studied through the two-particle correlations [45j 


cy k , = I d 3 ri<j 3 r 2 /(K, ri)/(K, r 2 )S (r - (ri - r 2 )) 
^ ' 1 f d 3 r 1 d 3 ?' 2 f(K,r 1 )f(K,r 2 ) 


(14) 


where /(K, r) is the particle phase space density in the final state. The information about the correlations is extracted 
through fitting a Gaussian form to the function S, 


S( K. 


r oc e 


(15) 


defining the femtoscopic radii i? ou t, -Rside, -Riong- The results for these extracted radii for pions are shown in Fig. [7] 
for d+Au, ,3 He+Au, p+Pb and Pb+Pb collisions, comparing hydrodynamic and non-interacting evolution. From 
this figure, one can observe a striking similarity for all the extracted radii between strongly interacting evolution 
(hydrodynamics) and non-interacting evolution (free streaming) for all simulated systems, small and large. Similarly 
to what was found for the case of radial flow, the femtoscopic radii are essentially insensitive to the details of the 
system evolution, as long as energy and momentum are conserved. 

In essence, this disqualifies the use of pion femtoscopic measurements as serving as evidence for a hydrodynamic 
phase during the system evolution. 
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FIG. 7. Pion femtoscopic radii (“-R ou t, Aside, -Ri ong ”) from simulations of granular d+Au, 3 He+Au, p+Pb and Pb+Pb collisions. 
Shown are results for identified particles for free-streaming evolution (no-interaction) and almost ideal hydrodynamics (r]/s = 
0.08), followed by a hadronic cascade. See text for details. 


D. Inverse slope parameters 


In the above result, it was found that radial flow is a feature of both almost ideal hydrodynamics and free-streaming. 
One might therefore be suspicious that the ubiquitous presence of radial flow in both models indicates a systematic 
failure of the modeling procedure since it is known from experimental data that radial flow does disappear in ’low’ 
multiplicity p+p and p+A collisions. 

Thus it is interesting to study if radial flow persists in the almost ideal hydrodynamic and free-streaming models 
if studying low-multiplicity p+A collisions. To this end, p+Pb collisions at yfs = 5 TeV were simulated for various 
multiplicity bins. For each multiplicity bin, the particle spectra for pions, kaons and protons were determined and fit 
with a form proportional to exp \Jp? + /T eS ^j with /x the pion, kaon and proton mass, respectively. The effective 

slope parameter T eS is reported in Fig[8^ along with the corresponding parameter measured for p+Pb collisions at 
yfs = 5 TeV by the CMS experiment [52] 7 Quantitative results can not be expected to match because the simulation 
results use an equation of state very different from that of QCD, but qualitative trends should be robust. The 
experimental data shown in Fig. [8] clearly shows that the effective slope parameter dependence on mass decreases 
significantly from the highest multiplicity selections to the lowest ones. This is often interpreted as a breakdown of 
collective behavior, as it can be linked to the disappearance of radial flow. 

In the simulation results, especially in the free-streaming model, one can recognize the same qualitative trend: the 
effective slope parameter dependence on mass decreases for lower multiplicity p+Pb events because the system does 
not ’live’ long enough to build up significant amounts of radial flow. Thus the simulation results are not inconsistent 
with the experimental findings. 

However, within the simple geometric Glauber model used for the initial conditions (see section IIA I, it is not 
possible to realistically describe either the highest multiplicity events (0-1% or higher) or the lowest multiplicity events 
(95%-100% or lower). In the Glauber model, it is not possible to have less than one collision, which according to 
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FIG. 8. Effective slope parameter from pion, kaon and proton spectra. Left and middle: Simulation results for almost 
ideal hydrodynamics and free-streaming for 0-10%, 20-30%, 50-60% and 90-100% multiplicity bins. Right: Experimental 
measurement by CMS experiment (figure from Ref. [52J) in terms of total particle tracks (8 to 235). 


section p I A| thus leads to a fixed amount of energy deposition, effectively putting a lower limit of the total multiplicity 
that can be simulated. Thus, the striking change seen in the experimental data from events with the highest and 
lowest number of particles cannot be simulated within the simple model for initial conditions adopted here. 


IV. SUMMARY AND CONCLUSIONS 

In this work flow signatures arising from two very different dynamics in the hot QCD phase following relativistic 
ion collisions very studied. In the first case, the hot phase dynamics was assumed to be described by non-interacting 
particles. In the second case, the hot phase dynamics was assumed to be described by extremely strongly interacting 
modes leading to almost ideal hydrodynamics. In both cases, the exact same initial conditions were implemented 
and the dynamics was required to correspond to the same equation of state. Also, in both cases the resulting 
energy-momentum tensor information was recorded on the same space-time grid and then passed on a hadron cascade 
“afterburner” using the same switching procedure. 

Because of this procedure, the resulting particle spectra between the two extreme cases of non-interacting and 
strongly interacting hot phase dynamics are directly comparable, and lead to the following findings: 

• Non-interacting (free-streaming) particle dynamics generally leads to equal or larger radial flow than strongly 
interacting dynamics (hydrodynamics) in all systems considered (d+Au, 3 He+Au, p+Pb, Pb+Pb). Also, as 
demonstrated in Fig. [2] radial flow velocities in hydrodynamics and non-interacting dynamics are similar already 
in the hot QCD phase, suggesting that this result is not dependent on the details of switching or the hadronic 
cascade. The overall amount of radial flow generated seems to be proportional to the time the systems spend in 
the hot QCD phase, naturally explaining why radial flow is observed to be very small in e.g. p+p collisions at 
y/s = 200 GeV (cf. Ref. [53]), which have a very short lifetime. This strongly suggests that the presence 
of radial flow extracted from experimental measurements should not be used as an indication for 
the presence of a hydrodynamic phase. 

• Non-interacting (free-streaming) particle dynamics generally leads to femtoscopic radii i? ou t, -Rside, Along that are 
very similar to those found in strongly interacting dynamics (hydrodynamics) in all systems considered (d+Au, 
3 He+Au, p+Pb, Pb+Pb). This strongly suggests that the results from femtoscopic measurements 
should not be used as an indication for the presence of a hydrodynamic phase. 

• Non-interacting (free-streaming) particle dynamics generally leads to considerably smaller elliptic flow than 
strongly interacting dynamics (hydrodynamics) in all systems considered (d+Au , 3 He+Au, p+Pb , Pb+Pb ). 

This strongly suggests that the presence of a sizable elliptic flow component extracted from 
experimental measurements is indicative of a hydrodynamic phase. 

• Non-interacting (free-streaming) particle dynamics generally leads to triangular and quadrupolar flow compo¬ 
nents that are comparable or even larger than hydrodynamics in light-on-heavy-ion collisions (d+Au , 3 He+Au and 
p+Pb). This suggests that higher order flow components extracted from experimental measurements for these 
small systems are not indicative of a hydrodynamic phase during the system evolution. It also suggests that the 
use of higher order flow components V 3 , tq, V 5 ,... as a high-precision “viscometer” in both small and large sys¬ 
tems should be reconsidered because non-hydrodynamic contributions can lead to a considerable contamination 
of extracted viscosity values. 
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As an outlook, one should note that generalizations of the free-streaming model description employed here, notably 
implementations of a QCD equation of state and weak, but non-vanishing interactions are possible should a direct 
comparison to experimental data become desirable. 
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Appendix A: Exact correspondence example between free-streaming and ideal hydrodynamics 


The similarities between non-interacting dynamics (free-streaming) and almost ideal hydrodynamics discussed in 
the main text may come as a surprise to some readers. However, both dynamics are just different formulations of 
energy-momentum conservation, so it may not be too surprising to find many similarities between these very different 
approaches. 

To elucidate the power that energy-momentum conservation places on the dynamics, let us give an example that is 
analytically solvable in both cases: that of SO(3) symmetric flow in Minkowski space (see Ref. [54] for more details). 

To wit, let us study equilibrium initial conditions at to = 0 for a system having a conformal equation of state 
Tff = e — 3 P = 0 with a spherically symmetric initial energy density profile given as 


e(t = 0, r) 


16L 4 

(A 2 + r 2 ) 4 ’ 


L = const. 


(Al) 


Using spherical coordinates r, 9, <\> one finds a metric tensor g ab = diag(l, — 1, — r 2 , — r 2 sin 2 0) and the associated 
non-vanishing Christoffel symbols as 

I 'ee = ^r, T^ =-r sin 2 6>, T 9 g = ^ , T^ = - cos 6 sin 9, T^ = ^ , T^ = cot6». (A2) 

In these coordinates, the equations for energy momentum conservation are V g T^ 1 ' = d^T^ + T^ la T al ' + = 0. 

Since the problem is spherically symmetric this implies T 4 ^ = T t9 = T r $ = T rS = 0, T 9 = and all components of 
T independent of </>. The equations of motion then simplify to 


dt (r 2 T 44 ) + d r (r 2 T rt ) = 0 , 
dt ( r 3 T rt ) + d r ( r 3 T rr ) - r 2 T 44 = 0 , 

d g T 99 = 0. 


(A3) 


Decomposing the energy-momentum tensor as in Eq. §, and putting n = 0 because of symmetry, one finds that a 
particular solution to the equations of motions and initial conditions is found as [54] 


e(t,r) = ~ 


t 2 L 2 + (LW*)’ 


u r (t,r) = 


rt 




7T^ = 0. 


(A4) 


It should be stressed that in this case the shear stress tensor turns out to be zero for all times because the shear 
velocity gradient is vanishing exactly, so that the present solution is exact for arbitrary values of viscosity. 

If one was to calculate spectra for massless particles with Boltzmann statistics at some switching hypersurface 
T = Tsw , the hydrodynamic solution would imply 


/hydro (t, T, p) 


Z exp 


—p°u° + p r u r 
T S W 


Z exp 


— |p| (A 2 +r 2 +t 2 — 2rf cosy) 
T S wVW r + (£2 +r 2_ 


(A5) 
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where p r = |p| cosx- The switching hypersurface in this case is located at those points (f,r) where e(t,r) = const. 
Choosing the proportionality constant so that e = T, we have 


'hydro 


(t) = t 2 - A 2 + 


2 , ^L\J\ — t 2 Tg W 


Tsw 


Tsw <2/A. 


(A6) 


Let us now solve the same problem with the same initial conditions for the non-interacting particle case. The initial 
particle distribution function is given by 


f F s{t = 0, r, p) 


Z exp 


~ —|p| (L 2 + r 2 Y 
2 A 


(A7) 


For Minkowski-space, the free-streaming solution to the Boltzmann equation is easily found to be f(t, r, p) = /(|x — ^ |), 
so that with above initial condition this leads to 


fFs(t,r, p) = Z exp 


— |p|(A 2 + r 2 + t 2 — 2rf cosx) 
2A 


(A8) 


With the choice for Z that is consistent with e = T 4 from above, one finds that the above solution corresponds to an 


energy density and flow velocity in Eq. (A41. Thus, for a switching hypersurface at T = Tsw one finds 


' FS 


(^) ^hydro 


{t) = t 2 — L 2 + 


2ZV1 - t 2 T 2 


sw 


Tsw 


Tsw < 2 /L. 


(A9) 


At this switching hypersurface, we thus have 2 L = Tsw \/Tt 2 L 2 + ( L 2 + r 2 — t 2 ) 2 and thus the free-streaming particle 
spectrum is 


fFs{t,r, p) = Z exp 


0 0 i TV 

-p u + p u 


Tsw 


= Z exp 


— |p| (A 2 + r 2 + t 2 — 2r t cos y) 
T sw V / 4t 2 A 2 + (A 2 + r 2 - 


/hydro • 


(A10) 


Thus, the free-streaming dynamics, including the particle spectrum at the hypersurface, is identical to the ideal 
hydrodynamic result. 

In particular, this proves analytically that free-streaming dynamics (no-interactions, no coupling to hadronic cas¬ 
cades) generates radial flow, as the particle spectra ( A10|A5 1 are equal and ideal hydrodynamics does generate radial 
flow. 


Appendix B: Importance of full stress-tensor matching for freeze-out 

As pointed out in the main text, it is important that the full energy-stress tensor is matched when switching from 
the hot phase to the hadron gas phase. An example of this can be given through monitoring the time-evolution of 
the momentum anisotropy, defined as 


/ d 2 x±r xx - Tyy 
p ~ f d 2 x±T xx + Tvv ’ 


(Bl) 


and its ’ideal hydrodynamic’ approximation obtained from dropping all viscous stresses from T a 


b. 


ideal _ 

e p — 


J d 2 x_ j_(e + P) ( u x u x — u v u v ) 

J d 2 x j_(e + P) {u x u x + uvuv) + 2P 


(B2) 


For the case of the smooth collision geometry considered in section |II E[ the time evolution of these quantities for 
free-streaming is plotted together in Fig. [9] with the evolution of the spatial anisotropy 


/ d 2 x±e ( x 2 — y 2 ) 
J d 2 x±e ( x 2 + y 2 ) 


(B3) 


From Fig. 
(as it should 


[blone finds that while the full momentum anisotropy e p is consistent with zero during the entire evolution 
be in free-streaming), this is not the case for the ideal hydrodynamic approximation ep deal . Thus, matching 



























17 


Anisotropy Evolution for Free-Streaming 



FIG. 9. Time evolution of spatial and momentum anisotropies in free-streaming for smooth initial geometry defined in Sec|II E| 


only the ideal hydrodynamic part of the energy-stress tensor in a freeze-out procedure would lead to a ’fake’ momentum 
anisotropy (fake elliptic flow). This has been observed in Ref. [47] , By contrast, the matching procedure in Eq. (12) 
has been designed to match the full energy-stress tensor, thus preventing fake contribution to arise from the hot phase 
to hadron gas switching procedure. 

Also, the example at hand demonstrates that the net zero momentum anisotropy in free-streaming is a consequence 
of a delicate cancellation between the (net positive) contribution in the ideal hydrodynamic part ep deal and a (net 
negative) contribution from the viscous stresses. 
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